Dissipative vs. Conditional Generation of Gaussian Entanglement and Spin Squeezing 
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Spin squeezing of collective atomic spins can be achieved conditionally via probing with light and subsequent 
homodyne detection, as is done in a Quantum Nondemolition measurement. Recently it has been shown that 
squeezing can also be created unconditionally by a properly designed dissipative dynamics. We compare the 
two approaches in a Gaussian description, and optimize over all Gaussian light-matter interactions. We find 
that in the optimal unconditional scheme based on dissipation the level of squeezing scales with optical depth 
as d -1 / 2 . In contrast, the optimal conditional scheme based on measurement of light - which in fact is not a 
Quantum Nondemolition measurement - can provide squeezing which scales as dT 1 . Our results apply directly 
also to the creation of entanglement in the form of non-local spin squeezing of two atomic ensembles. 



I. INTRODUCTION 

Engineered dissipative dynamics can stably generate non- 
classical states and nontrivial dynamics of quantum systems, 
as was shown in recent theoretical [1-17] and experimen- 
tal [18, 19] work. Dissipative generation of spin squeezing 
[1], entangled states of single atoms [2-4], atomic ensem- 
bles [5-10, 18] or micromechanical oscillators [11, 20], and 
quantum states of many body systems [12, 13] is possible, 
and even quantum computation [14], quantum communica- 
tion [15], and quantum simulations [16, 19] can be performed 
purely based on dissipation. Another method for engineering 
of quantum states and dynamics, which likewise rests upon 
the coupling between the system to be controlled and its en- 
vironment, is continuous measurement and feedback control 
[21-26]. Here the state or dynamics of the system is condi- 
tioned upon a measurement performed on the environment, 
and may be controlled via feedback depending on the mea- 
surement outcomes. Continuous measurement and feedback 
control is a well established method in an equally broad range 
of physical systems [18, 27-39]. This naturally raises the 
question of which methods work best for a given task or figure 
of merit. 

In this article we address this question for the case of 
spin squeezing [40] and spin squeezing based entanglement 
of atomic ensembles [41]. The creation of spin squeezing 
can be considered as the prime example for measurement 
and feedback based quantum control. It relies on the optical 
probe of a collective atomic spin by off-resonant light, homo- 
dyne detection of light, and an appropriate feedback on the 
atomic spin [42]. The light-matter interaction is thereby com- 
monly tailored such that a quantum nondemolition measure- 
ment (QND) [43, 44] of the atomic spin is realized. Squeez- 
ing of atomic spins via (continuous or pulsed) QND probes 
with and without feedback was demonstrated in [42, 45-51]. 
When squeezing of two collective spins is realized in this way, 
entanglement in the form of Einstein-Podolsky-Rosen (EPR) 
squeezing is established among the two atomic ensembles 
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FIG. 1. (Color online) One (a) or two (b) bosonic modes (collective 
atomic spins) couple sequentially to ID field, which might be subject 
to homodyne detection after the interactions. 



[41], as was done in [52]. In the alternative context of dis- 
sipative quantum state engineering it was shown recently that 
a careful choice of the light-matter interaction can provide a 
spin squeezed or entangled state in steady state without apply- 
ing any measurement on the light and feedback on the atoms 
[7, 8]. This was realized experimentally soon after [18], albeit 
with the help of moderate feedback compensating for experi- 
mental imperfections. 

Prompted by these theoretical studies and experiments we 
investigate here how the light-matter interaction has to be tai- 
lored in order to achieve a maximal amount of spin squeez- 
ing or entanglement. We consider and compare all three rel- 
evant scenarios: unconditional, dissipative state generation, 
measurement based conditional state generation, and uncondi- 
tional measurement and feedback assisted state generation. In 
each case we optimize over the class of Gaussian light-matter 
interactions (i.e. those which are quadratic in (collective) cre- 
ation and annihilation operators, see below). Crucially, we 
take into account decoherence induced by spontaneous emis- 
sion due to the off -resonant light probe. We do this on the ba- 
sis of a minimal, yet microscopically justified Gaussian model 
[53]. This approach allows us to relate the spin noise reduc- 
tion or entanglement achievable under optimized conditions 
eventually to the optical depth d of the atomic ensemble [54]. 

Our results are as follows: First, for unconditional, dissi- 
pative state generation the optimal (EPR) squeezing scales as 
dr 1 / 2 . This is the same scaling as is achievable with a QND 
measurement [55]. While the scaling with optical depth is 
the same in both cases a QND measurement provides squeez- 
ing which is twice as large for the same parameters. How- 
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ever, QND squeezing is prepared conditionally and becomes 
unconditional only when feedback is applied. Thus, dissipa- 
tive generation of squeezing spares the need for feedback at 
the cost of a 50% reduced level of squeezing. Second, mea- 
surement induced squeezing can exhibit a rather remarkable 
scaling of dr 1 for the regime of moderate optical depths. In 
this most relevant regime we thus predict an enhancement of 
dr 1 ! 2 in comparison to the scaling for QND squeezing. For 
large optical depths the optimal scaling eventually approaches 
dr 1 / 2 again. The enhancement comes at the cost of the ne- 
cessity of applying feedback in order to stabilize the spin state 
as the conditional state is unstable (in contrast to the case of 
QND squeezing). In each case we determine the optimal light 
matter interaction as a sum of beam splitter and two mode 
squeezing dynamics, with dominating beam splitter contribu- 
tion for the dissipative generation, and dominating two mode 
squeezing for the measurement based scheme. 

Our results are formally derived by means of Gaussian 
quantum stochastic master equations and corresponding equa- 
tions for displacement vectors and covariance matrices {i.e. 
the first and second moments) [56, 57]. We supply the basics 
and some useful relations for this formalism in Appendix A. 

The paper is organised as follows: In Sec. II we introduce 
a model describing atomic ensembles interacting with an off- 
resonant probe light. In the next Section III we derive the 
stochastic master equation for a single bosonic mode probed 
by light. After that the corresponding equations for variances 
of the oscillator quadratures are introduced and solved for 
conditional and unconditional dynamics. In Sec. IV we con- 
sider two bosonic modes (spin ensembles) coupled to a ID 
field in a cascaded fashion. We show that the evolution of 
EPR modes of the two oscillators coincides with the dynam- 
ics of the single bosonic mode considered in previous section. 
The discussion of our main results for spin squeezing and en- 
tanglement is given in Section V. 



II. MODEL 

We want to describe a general interaction of an atomic en- 
sembles of spin polarised atoms with a ID light field which 
leads to spin squeezing or generation of entanglement. The 
widely used procedure for a macroscopically big collective 
spin is to approximate the spin component along its axis of 
polarisation with its mean value J x m (J x ). Rescaling of the 
rest of the spin components provides the canonical position 
X = Jy/y {Jx) an d momentum P = J z /y/ (Jx) operators 
with a proper computational relation [X, P] = i. The rig- 
orous description of this procedure is given by the Holstein- 
Primakoff transformation [58]. It allows us to work with 
bosonic modes and obtain results which can be applied to 
macroscopic ensembles of spin polarised atoms. 

The optical probe will inevitably cause decoherence in the 
spin systems, which can be taken into account by introducing 
decay rates and effective occupation numbers of the environ- 



ment: 

X = - lx X + f x (t), (f x (t)f x (t')) = n x S(t - t 1 ), 

P = -1 P P + f P (t), (f P (t)f P (t')) = n p S(t - t'). 

The relevant decay rates and noise terms due to the interaction 
of atoms with the reservoir of electromagnetic modes can be 
calculated form first principles, as was done in [53]. Deco- 
herence from other sources (if present) can be added indepen- 
dently. In order to keep the following discussion conceptually 
simple we will assume equal decay rates ^ x = j y = 7 and oc- 
cupation numbers n x = n y = n. Our approach can be easily 
adapted in order to account for more general situations. The 
equivalent description of the decoherence process in terms of 
the master equation is given by 

p = 1 {n + l)V[a]p + 1 nD[a)]p, (1) 

where the Lindblad term is defined as T>[x]p = xpx^ — 
\{x^x,p} + . 

We consider a setup as shown in Fig. 1. A ID field is 
coupled to one or two localized bosonic modes positioned at 
points z\ and Z2- In the latter case the field thus interacts se- 
quentially with the two localized modes at times ti = zi/c 
and t2 — z-ijc > t\. After the interaction the field can be ho- 
modyne detected. We are interested in the stationary state of 
the localized bosonic modes, and in particular in its squeezing 
and entanglement properties. Stationarity is thereby under- 
stood with respect either to the purely dissipative, uncondi- 
tional dynamics or to the dynamics conditioned on the homo- 
dyne detection. 

III. A SINGLE OSCILLATOR COUPLED TO A ID FIELD 

Conditional Master Equation 

We consider first a system with one localised bosonic mode 
coupled to a ID field as shown in Fig. la. The ID field is 
described by bosonic operators in frequency space b^, or in 
time 




where t = z/c. The corresponding quadrature operators are 
x{t) = [b(t) + 6(t)t]/V2 and p(t) = -i[b(t) - b{tf]/V2, 
and carry dimension Hz 1 / 2 . The oscillator is described by a 
dimensionless bosonic operator a with corresponding canoni- 
cal operators X = (a + a^)/\/2 and P = —i(a — aJ)/y/2. 

We assume a general quadratic interaction of the field with 
the oscillator. Without loss of generality [59] any such inter- 
action can be written as 

H mt = ^[aXx(t)+(3Pp(t)], 

where g has dimension Hz and determines the strength of the 
coupling, and the coefficients a and f3 can be parameterized 

as 

a = cos(6>) P = sin(0) (2) 
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The ratio of these coefficients determines the qualitative char- 
acter of the coupling. When written in terms of creation and 
annihilation operators the Hamiltonian becomes 



where 



V2 



(aX + i/3P) 



a + (3 a- (3 



Thus, for 9 = the Hamiltonian describes a QND interaction 
~ Xx(t), while for 9 = tt/A it is a beam splitter interaction^ 
atf (t) + aJb(t), and for 9 = —n/4 it is a two mode squeezing 
interaction ~ ab(t) + a'b'(t). Other values of 9 cover any 
quadratic interaction in between. 

The full Hamiltonian for the oscillator and the ID field is 



H = H + H] 



bath 



H" 



where H is a local Hamiltonian for the localized mode, and 
-Hbath = /dww&^b^. Eliminating the bath in a Born-Markov 
approximation yields a master equation 



p = -i[H,p]+gV[s]p (3) 

where the Lindblad term is defined as T>[x]p = xpx* — 
\{x i x 1 p}+. 

To the master equation (3) we can now add other processes 
relevant to the dynamics, as explained in Sec. II: First, we 
assume that the oscillator couples at a rate 7 to a thermal 
reservoir, causing in thermal equilibrium a mean occupation 
number of n, and therefore implies a thermal decoherence rate 
ryth _ Second, jf the transmitted field is subject to homo- 
dyne detection of both conjugate light quadratures (after being 
split on a beamsplitter with reflectivity e, cf. Fig. 1) the con- 
ditional state evolves according to the conditional stochastic 
master equation 



dp 



- i[H,p]dt 



gV[s]pdt 
TjHlse^pdWi 
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r giH[ise^]pdW 3 
l)V[a]p + jnV[a^]p} dt, (4) 



where H[x]p — (x— (x))p+p(x> — and dWi is a Wiener 
increment of zero mean and dW 2 = dt. <fr is the phase of 
the local oscillator in the homodyne detection. By convention 
= corresponds to a measurement of the phase quadrature 
Pit)- 



Steady State Variances 

For Gaussian states the stochastic master equation (4) im- 
plies a time evolution which is described by a stochastic linear 
differential equation for the conditional mean values (X(t)) 
and (P(t)), and a deterministic nonlinear differential equation 
for the conditional variances 



V c = AX< 



U c = AP A 



(5) 



and covariances C c = (XP + PX) - 2(X) (P). In [17, 56] 
compact formulas were derived in order to pass from a general 
stochastic master equation such as (4) to the corresponding 
equations for V c , U c , and C c . In Appendix A we provide a 
summary of this extremely useful formalism. 

For an optimal choice <p = of the local oscillator phase 
(giving rise to the strongest squeezing) and a beamsplit- 
ter reflectivity e the conditionally squeezed variance V c and 
the conditionally unsqueezed variance U c obey the nonlinear 
equations of motion 

V c = - [7 - .9(1 - 2e)ap]V c - g(l - e)a 2 V c 2 

+ 7 (2n + l)+ 3 e/3 2 , (6) 

Uc = - [7 + 5(1 - 2e)aP]U c - geP 2 U 2 c 

+ 7 (2n + l)+ 5 (l-e)a 2 , (7) 

as follows from Eqs. (A4) in Appendix A. We suppress the 
corresponding equation of motion for the conditional covari- 
ance C c , as it will vanish in steady state for the present case. 
The conditional variances have the steady state solutions 



Vr. 



1 



2(1 - e)ga 2 



{(l-2e) ff a/?- 7 



Uc = 



+ \/irt - 9 a PY + *-yga[(2n + 1)(1 - e)a + e/3]}, (8) 
l —{-{l-2e)gaP- 1 



2eg(3i 



+ V(7 + 9(*P) 2 + 4jg(3e[(2n + l)/3 - a] }. 



(9) 



The result for the conditional variances in steady state (8) 
and (9) can be compared to the unconditional variances V u 
and U u , which can be obtained by solving (4) neglecting the 
stochastic term proportional to the Wiener increments. Ap- 
plying again formulas (A4) the master equation implies for 
the unconditional squeezed and unsqueezed variances 



K = -(7 + 9<*P)V U + gP 2 + i(2n + 1), 
U u = -(7 + gaf3)V u + ga 2 + 7 (2n + 1) 

with the steady state solutions 



U u 



,g/3 2 + 7 (2n + 1) 
7 + ga/3 

ga 2 +7(2n+ 1) 
7 + 9af3 



(10) 

(11) 

(12) 
(13) 



Equations (8), (9), and ( 12) are the main results of this section. 

Instead of asking for the steady state corresponding to a 
continuous wave drive and a continuous detection of light, we 
can also ask for the solution to the equation of motion for 
the conditional variance V c for finite times corresponding to 
driving the system with a pulse. Integrating (6) for a time r, 
measuring only one quadrature of light (e = 0), neglecting any 
losses (7 = 0), and selecting a QND interaction (9 = 0) with 
initial condition V c (0) = 1 yields V c (r) = 1/(1+<?t). We can 
compare this to the known result for QND squeezing, which 
is given by V^9 nd (t) = 1/(1 + n 2 ) where k 2 = dg with 
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optical depth d and 77 = 7T the degree of atomic depumping 
from a pulse with duration r, see e.g. [41]. We can therefore 
conclude that 



.9 = <fy 



(14) 



Since the expressions for conditional (8), (9) and uncondi- 
tional (12), (13) variances in steady state depend only on the 
ratio 7/17 the results of the substitution g = dj become in- 
dependent of the decay rate, and depend only on the optical 
depth, the effective mean thermal occupation, and the param- 
eters e and 9 characterzing the measurement and the light mat- 
ter interaction. Substituting (14) and (2) into (8), (9) and (12), 
(13) we arrive at the final results given in Sec. V, Eqs. (18) and 
(19). Before we enter into the discussion of these results we 
will show that they apply directly also to the EPR squeezing 
of two bosonic modes. 



IV. CASCADED COUNTERROTATING CONFIGURATION 

The preparation of a squeezed state for an oscillator is 
closely related to the production of an entangled state of two 
oscillators. Essentially one needs to prepare EPR modes of 
the oscillators in a squeezed state. In order to achieve this we 
consider two oscillators coupled to the same field mode in a 
cascaded fashion as shown in Fig. lb. The full Hamiltonian 
for the cascaded system of the two localized modes and the 
ID field is 



H — H 1 + H 2 + -ffbath + H[ n 



m 



where Hi are local Hamiltonians for the two modes, and 

= ^[atXixiti) + ftPiP(ii)] , 

with t 2 > t\. To ease notation we assume couplings of equal 
strength <?j=i,2 = g- Eliminating the bath in a Born-Markov 
approximation yields a cascaded systems master equation [4, 
60-63] 



P 



Hi + H 2 - «|(4 s i _ 4 S 2), P + g1>[si + s 2 ]p 

(15) 



where the Lindblad term depends on a collective jump term 
s\ + s-2- The effective interaction among the two modes me- 
diated via the ID field is 



•9 1 t f \ 
-i-{s 2 Si - s[s 2 ) 



■fj 

l- 
2 



a\a 2 — a[a 2 



a 2 /3i + aij3 2 ( + + 
I aia 2 - a[a 2 



Note that if both modes couple in the same way to the ID 
field, a\ — a 2 , f3i — f3 2 , this mediated interaction is a pure 
beam splitter like coupling. 

As we did above with a single oscillator, we can now add 
to the cascaded systems master equation (15) other processes 
relevant to the dynamics: First, we assume both modes couple 



at a rate 7$ to a thermal reservoir causing in thermal equilib- 
rium a mean occupation number of rii, and therefore implying 
a thermal decoherence rate 7*" = r fin i . Second, if the trans- 
mitted field is subject to homodyne detection the conditional 
state evolves according to the conditional stochastic master 
equation 



dp 



Hi + H 2 - i|(4 s i - 4 S 2), P 



dt 



+ gV[si + s 2 ]pdt + y/g(l-e)H[( Sl + s 2 )e^]pdWi 

+ y /gZH[i(a 1 + 8 2 )e i *]p&W2 

+ {rt n i + + iniV[4]p} dt. 

i 

We now specialize this fairly general model. First, we as- 
sume the "Copenhagen setup" of counterrotating harmonic 
oscillators in which a broad range of quantum information 
protocols were performed [41]. In this setup the local Hamil- 
tonians are 

H\ + H 2 = wajfli — Loa\a 2 . 

such that one of the ensembles effectively realizes an oscil- 
lator of negative mass. In a rotating frame with respect to 
Hi + H 2 one finds 



Sl (i) + a 2 (t) =s + e~ iu t + s^ t 



where 



oti + Pi 
s+ = — 2— a x 

a 2 + [3 2 
s- = r a 2 



a 2 - (3 2 + 

n 1 



oil - Pi i 
r a . 



Assuming g, jf 1 <C w we can perform a rotating wave ap- 
proximation such that the conditional master equation be- 
comes 



dp 



-ig{a 2 fii - aif3 2 ) 



(aia 2 - a\a\),p 



dt 



-T>[s+ + s-]pdt + 1 2? [s_|_ — S-]pdt 



H[{s + + S-)e^}pdW c i 




U[i{s+ + s.)e^]pdW c2 

U[-i{a + -sJ)e i +]pdW, x 
U[{s+ - s.)e^]pdW s2 
+ {^ ni + l ) v \. a i]P + Hn&WAp} d< - 

i 

Here dW / c i( 2 ) and dW^m) are independent Wiener incre- 
ments, corresponding to noise in sine and cosine modulation 
modes at the sideband frequencies ±w for the x (p) quadra- 
ture. For details we refer to [57] (Sec. 4.5). 
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For an entirely symmetric situation a* = a, fa — fa 7j = 
7, and rii — n for both i = 1 , 2, we get 

s + + s_ = al + + i/3P + 
s + — s_ = j(a_P_ — ifaXJ) 

where the EPR operators are Jf± = (Xi ± X 2 )/V2 and 
P± = (Pi ± P 2 )/V2. In this case the dynamics factorizes 
in these modes, and the reduced density operators p± for the 
EPR modes fulfill the stochastic master equations 



TMS 



d P+ = ^V[aX + + ifaP + ] P+ dt 



(16) 



\/ g{1 2 e) n[(ax + + ifap + y% + dw cl 



-H[i{aX+ + ifaP+)e^)p+dW c2 

+ 7(n + l)V[a + ]p + dt + in V[a) + }p + dt, 
dp_ = ~V[aP- - ifaX_]p_dt 



(17) 



.9(l-e) 




«[-t(ctP_ - ifaX^)e l *]p^dW sl 



y«[(aP_ - i J 8A-_)e , «]p_dW, 2 
+ 7(n + l)D[a_]p_dt + 7nZ>[aL]p_dt, 



where a± = (X± + iP±)/v2. One can see that the equa- 
tions (16) and (17) obtained for the two EPR modes (+ and 
— ) are equivalent to the stochastic master equation for a sin- 
gle oscillator considered in Sec. (Ill), cf. (4). Therefore, the 
steady state solutions for the single oscillator's variances (8), 
(9), (12) describe the squeezing of EPR modes as well. En- 
tanglement of the two oscillators is established if the sum of 
the two squeezed EPR variances falls below two units of spin 
projection noise (that is, 2 x 1/2 in our convention). 



V. DISCUSSION 

Now we come to the discussion of the results obtained for 
the squeezing of a collective spin or entanglement of two of 
them. As explained at the end of Sec. Ill the expressions for 
the unconditional (12) and conditional (8) variances in terms 
of the interaction parameter 9 and the optical depth d read 



2n + 1 + ds'm 2 9 
l + rfsin(26»)/2 ' 



(18) 



Vc. = 



1 



2(1 -e) 



(1 + e) tan0- 



1 



dcOS 2 9 




1 



dcos 2 9 



tan 9 



4[(2n + l)(l 



e tan ( 



dcos 2 9 



(19) 
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FIG. 2. (Color online) Unconditional squeezed variance V u 
(Eq. (18)) versus interaction parameter 6. The dynamics is stable 
only for 8 > 8 C where 6 C is given in (20). V u exhibits a minimum 
for an optimal choice of #° pt given in (21) in a region where the 
beam splitter (BS) interaction dominates over the two mode squeez- 
ing (TMS) interaction. Curves for the following parameter values are 
shown: d — 5, n — is the blue (dark grey) line, d — 50, n — is 
the green (light grey) line, and d — 50, n = n c is the red (top) line. 



Unconditional variance 

• The unconditional steady state exists only provided 9 > 
9 C where the critical interaction parameter is given by 

7-r-#sin(20 c )/2 = 0, or 



= — - arcsin(2/d). 



(20) 



Otherwise the system becomes dynamically unstable, 
see (10). However, this is not the case for the condi- 
tional state which can be stabilised for any interaction 
parameter 9 by the measurement of both quadratures of 
the output light (e ^ or 1) and a feedback (see below). 

• The unconditional variance (18) exhibits a tradeoff with 
respect to 9, as is shown in Fig. 2. The minimum of 
the variance (maximum of entanglement) is achieved 
for the optimal interaction parameter 9° pt given in the 
Appendix B, and in general requires an interaction dom- 
inated by the beam splitter dynamics. In the limit of 
high optical depth it reads 



l9°pt 



2n+ 1 



(21) 



The exact value of 9° pt is shown in Fig. 3 as a function 
of optical depth. 

• Substituting the optimal interaction parameter in 
Eq. (18) one obtains the minimal unconditional vari- 
ance 



yopt = 2 



d^/(2n + l)(2n + 1 + d) + 1 -d- 2(2n + 1) 



d 2 -4 



2n + l 
d 



(22) 
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FIG. 3. (Color online) Optimal interaction parameter 9 versus optical 
depth. Solid curves show conditional (the lower blue line) 9° pt and 
unconditional (the upper green line) 8° pt optimal interaction param- 
eters. Dashed lines are asymptotics. For large optical depths the op- 
timal interaction for dissipative generation of squeezing approaches 
the QND interaction (9 = 0), while the best choice for conditional 
generation of squeezing becomes the two mode squeezing interaction 
(9 = -tt/4). 

where the approximation holds for d ^> 1. Fig. 7 be- 
low shows V° pt versus optical depth. We recover the 
scaling of squeezing as dr 1 / 2 which also applies for 
a QND measurement. This can be seen directly from 
(19) by setting 6 = and e = 0. In the asymp- 
totic regime for large d the QND squeezed variance is 
then V C QND = y/(2n + l)/d, half of the dissipatively 
squeezed one. 

Conditional variance (e = 0) 

• The conditional variance (19) exhibits a similar trade- 
off. It has a minimum with respect to 6, as shown in 
Fig. 4. The corresponding optimal interaction parame- 
ter 8° pt is defined in Appendix B and shown in Fig. 3, 
and in general corresponds to an interaction which is 
dominated by the two mode squeezing dynamics. It 
does not have a nice analytical form but in the limit of 
high optical depth can be expanded as 

0° pt ^-j + 2r ^+O(d- 2 ). (23) 

• From the expression (23) one can see that the optimal 
interaction for the maximum entanglement is close to 
the two-mode squeezing regime. The expression for the 
optimised conditional variance in the limit of d ^> 1 
reads 

Kopt ^2(2n + l) (24) 

The striking feature of the conditional variance (24) is 
the inverse scaling with the optical depth in contrast to 
the ~ 1/Vd scaling of the unconditional (22) squeez- 
ing (entanglement) for dissipative or QND generation 
for squeezing. The price we pay for this significantly 
improved scaling is the unbounded variance of the un- 
squeezed component. In the limit of e — > the ex- 
pression for the conditional variance U c (9) converges 
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FIG. 4. (Color online) Conditional squeezed variance V c (solid) 
and unsqueezed variance U c (dashed) versus interaction parameter 
6. The beamsplitter reflectivity is e = 0. Curves for the following 
parameter values are shown: d = 5, n — is the blue (dark gray) 
solid line for V c and dashed for U c , d = 50, n = is the green (light 
grey) solid line for V c and dashed for U c , and d = 50, n = n c is the 
red (top) solid line for Ve- 
to the unconditional unsqueezed variance U u (13) and 
therefore too becomes unstable for 9 < 9 C . 

Stabilised conditional variance (e > 0) 

• In order to profit from the enhanced scaling of squeez- 
ing with the optical depth the antisqueezed spin compo- 
nent therefore has to be actively stabilized via feedback. 
To do so it is necessary to detect also the second, con- 
jugate light quadrature, from which the required infor- 
mation about the antisqueezed spin component can be 
gained. This can be achieved in the measurement setup 
shown in Fig. 1. For a reflectivity e > a little bit of 
output light is substracted and subject to a homodyne 
detection of the amplitude quadrature. With appropri- 
ate feedback of the corresponding photocurrent the an- 
tisqueezed quadrature can be stabilized to a value which 
scales as U c ~ 1/e (in units of shot or projection noise), 
at the cost of a somewhat reduced squeezing of the other 
spin component. The variances for a small reflectivity e 
are show in the Fig. 5 and Fig. 6. 

• Most importantly, in Fig. 7 we see that the 1/d scaling 
of squeezing can be stabilized for optical depths up to 
a value d* ~ (2n + l)/e. For greater optical depths 
the scaling of squeezing is changed back to the inverse 
square root law V° pt {d > d*) ~ 2^(2n+ l)e/d. The 
important conclusion is that for a given optical depth 
d one can achieve an enhanced scaling of squeezing 
as 1/d by using a reflectivity e ~ (2n + l)/d. The 
corresponding level of antisqueezing will be given by 
U c (dJ) ~ \/2/e oc d. 

• In the regime of unconditionally stable dynamics 
(shown by a grey area in Fig. 5 and Fig. 6) the con- 
ditional state is dynamically stable too. Feedback may 
be applied in order to convert the unconditional squeez- 
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FIG. 5. (Color online) Squeezed unconditional variance V u (green 
(light grey)), conditional variance V c (blue (dark grey)) and un- 
squeezed conditional variance U c (blue (dark grey) dashed) versus 
interaction parameter 8 for an optical depth d — 5, minimal atomic 
noise n — and beamsplitter reflectivity e = 0.05. Interaction pa- 
rameter 9 — corresponds to a QND interaction, for < 9 < ~ the 
beam splitter interaction dominates, for — ? < 9 < the two mode 
squeezing interaction is dominant. The greyed region corresponds to 
the unconditionally stable dynamics. 
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FIG. 6. (Color online) Same as for Fig. 5 but for critical (effective) 
temperature corresponding to decay to a thermal state with mean oc- 
cupation number n = n c = d(v2 — l)/4. It is impossible to uncon- 
ditionally achieve spin squeezing or entanglement for this or higher 
occupation number. 



FIG. 7. (Color online) Squeezed and unsqueezed variances vs opti- 
cal depth for minimal noise n = 0: (blue (dark grey)) Conditional 
variance for the optimal 9° pt , (green (lower light grey)) Uncondi- 
tional squeezed variance for optimised interaction parameter 9° pt , 
(yellow (upper light grey)) Unconditional unsqueezed variance for 
optimised interaction parameter (9° pt , (red (upper)) Conditional un- 
squeezed variance for the optimal 9° pt . Dashed lines: (blue (dark 
grey)) Asymptotic scaling 2/d for conditional dynamics and (green 
(light grey)) 2/yd for the unconditional dynamics, (red (upper)) 
asymptotic of the conditional unsqueezed variance 1/e. 



of n which should not be surpassed in order to achieve an 
entangled state. 

• The unconditional steady state squeezing (V u < 1) can 
only be observed if the occupation number of the envi- 
ronment is below a critical limit n c = d{\/2 — l)/4. 

• The conditional state can remain squeezed (entangled) 
until the temperature reach n c ~ ^ — A + 0{dr z ) 
which is approximately six times higher than the uncon- 
ditional n c . 

• Yet another critical temperature for the conditional 
squeezing in the dynamically stable regime is given by 

0(d 
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CONCLUSION 



ing into conditional squeezing. Depending on the ap- 
plication it might not be necessary to do so, and the 
conditionally squeezed state might give the same result 
{e.g. in some Quantum Information protocol). On the 
contrary, the unstable regimes where the 1/d scaling of 
squeezing can be expected can be accessed only by ap- 
plying a continuous feedback. The master equation de- 
scribing the unconditional state of the system with feed- 
back is given in Appendix C. 

Finite occupation number If the decay does not go back 
to the ground state (fully polarized state) we can identify an 
effective temperature giving rise to a mean thermal occupation 
n. Fig. 6 shows the effect on the unconditional and conditional 
variances. For a given optical depth d there is a critical value 



We have considered spin squeezing and entanglement of 
collective atomic spins obtained in a dissipative dynamics. We 
performed the optimisation over all possible quadratic light- 
matter interactions in a Gaussian description. The optimal un- 
conditional squeezing exhibits a scaling with optical depth as 
dr 1 / 2 . We found that if one performs homodyne measure- 
ments of the probe light it is possible to achieve better squeez- 
ing with a remarkable scaling as d _1 by applying a feedback. 

Our results are based on a Gaussian description of the light 
matter interaction, and neglect the spatial extension of the 
atomic ensembles along the propagation direction. This is 
an excellent approximation for a QND interaction, but be- 
comes less appropriate for dominant beam splitter or two 
mode squeezing interactions. However, we expect that our ap- 
proximation still applies for ensembles at room temperature 
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where thermal motion averages out spatial inhomogeneities. 
In view of the possible gain in squeezing a more careful treat- 
ment based on integration of Maxwell-Bloch equations is cer- 
tainly desirable. 
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Appendix A: Stochastic Master Equations and Covariance 
Matrix Equations 

Consider a system composed of a bosonic modes which 
obeys a stochastic Master equation 



A P = -i[H, P ]At + ^V[Ji\pAt + ^H[J l ]dW l , (Al) 

i=l i=l 

where V[x]p = xpx* — \{x^x, p} + andH[x}p = (x—(x))p+ 
p(x* — (x*)). For a 2a vector of canonical operators r T = 
(xx,Pi, ■ ■ ■ ,x a ,Pa) me Hamiltonian 

H = - f T Mr, 
2 

is characterized by a real symmetric 2a x 2a matrix M, and 
the b jump operators 

are determined by b vectors ji £ C 2a . We assume that c out of 
the 6 decay channels are monitored, as described by the mea- 
surement terms proportional to the Wiener increments AWi. 
We restrict attention to the case where AWi are independent 
Wiener processes dW^dW^ = SijAt. 

If the state of the system is Gaussian it is fully described by 
the vector of first moments collected in the real 2a-component 
displacement vector s, and the real symmetric 2a x 2a covari- 
ance matrix T defined by 

Si = tr {pri} , Tij = tr {pinrj + r^)} - 2s l s j , (A2) 

such that for a single system and in the notation used in the 
main text 



r = 



(2V C C c 
IC C 2U C 



It is then possible to show that for Gaussian states the equa- 
tions of motion for s and T implied by (Al) are 

As = QsAt + (TA— aB)AW, 

f = (Q + 2<jBA t ^ T + T (Q T + 2AB T a T ^j 

P - 2aBB T a T ^j - 2rX4 T r. (A4) 

The real matrices which enter here are determined from M 
and ji as follows: 



Q = a(M + R), R = --J2(j*l ir -b.c.) 

i=i 

1 6 

P = 2aSa T , S = -J2(f:i T + h.c.) 

i=l 

and finally, a is the 2a x 2a symplectic matrix 

1 



' 



This can be proven by substituting (Al) into (A2), using the 
cyclicity of the trace, the canonical commutator [7»,rj] = 
icrij, and the property of Gaussian states that 

tr {p{r h {r m ,r n } + } + } = 

2 (F; m 5 n -\- T n ;5 m -(- T mn 5; -|- 2s n S n S m ) . 



Appendix B: Optimal 9 

An optimisation of the expression for the unconditional 
squeezing (18) with respect to the interaction parameter pro- 
vides the optimal 6>° pt given by 



(9° pt = arctan 



V(2ra + l)(2rt+l + d) + l- l 
2n + 1 + d 



2n+ 1 



(Bl) 



(A3) 



The approximation assumes d^S> 1. 

The conditional variance (19) optimisation boils down to 
finding a proper root of a cubic equation: 

= 2x 3 (l - 2e)(2n + l + de) 

- x 2 {4(2n + l) 2 - e[16n(n + 1) + d 2 } 
+ d(2n + l)[2-(l-2e) 2 ]} 

- 2x(l - 2e){2n + 1 - de) 
+ rf(2n + l)(l-2e). 

The optimal interaction parameter is then given by 6° pt — 
arctan(a;). 
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Appendix C: Feedback Master Equation 

Based on the light measurement outcome one can apply a 
Hamiltonian feedback to the oscillator or to the EPR modes of 
a pair of them. The feedback Hamiltonian is proportional to 
the measured photocurrents I\ (t) and I 2 (t) of the homodyne 
detectors, 



Hfb — Ii(t)Fi, 



(CI) 



where the feedback operators are F\ = £iP and F 2 = £, 2 X 
with £1 and £2 being the feedback gains. The resulting master 
equation reads 



- 2 [WaO- - e ) p i + Ws^^ 8 + h.c, P ] 

+ vy g {l - ejs - iF x ]p + geV[^gis - F 2 ]p, (C2) 



By choosing appropriate feedback gains £1 and £2 one can 
achieve steady state of the system with variances given by the 
conditional dynamics (8) and (9) for any interaction parame- 
ter 9. Applying Eqs. (A4) to the master equation (C2) we ar- 
rive at the Lyapunov equations for the variances Vfb and [/ft,: 



V fb = - 



j + ga/3 - 2a^2g(l - e) 
7(271 + 1) + gf - 2/?6 v/2g(l - e) + 2^, (C3) 



f/fb 



1 + ga0-2^2ge 



C/fb 



+ 7 (2n + 1) + ga 2 - 2a£ 2 ^2gt + 2£f. 



(C4) 
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